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I. INTRODUCTION 



The disparity between the energy scales of typical QCD phenomena and the scale for 
nuclear binding makes nuclear systems an ideal playing ground for Effective Field Theory 
(EFT) methods. The first motivation behind the EFT program in nuclear physics is to 
obtain a description of nuclear phenomena firmly based on QCD. Recently, considerable 
effort has been put into the development of a phenomenologically successful EFT for nuclear 
physics 1^, ^ ^. It soon became apparent that there is one feature which distinguishes the 
nuclear EFT from most of the many other applications of EFT methods: The low-energy 
expansion is non-perturbative in the sense that an infinite number of diagrams contribute 
at each order. This follows from the very existence of shallow real and virtual bound states 
whose binding energies are so small that they lie within the range of validity of the EFT 
even though they cannot be described by perturbing the free theory. In fact, there are 
at least in the two and three-nucleon sector states so loosely bound that their scales do 
not seem to be connected even to the soft QCD scales or Instead, they appear to 
require an additional "accidental" fine-tuning whose microscopic origin is not understood. 
However, this phenomenon offers the opportunity for a more radical use of EFT in the 
few nucleon case. In this approach, all particles but the nucleon themselves are considered 
high energy degrees of freedom and are consequently "integrated out". The resulting EFT 
is considerably simpler than potential models or the "pionful" version of nuclear EFT (in 
which pions are kept as explicit degrees of freedom), but its range of validity is reduced to 
typical momenta below the pion mass. Albeit this might seem a severe restriction, there 
are many processes situated in this range which are both interesting in their own right and 
important for astrophysical applications. Also, we believe that an understanding of this 
simpler EFT is a necessary condition for understanding more general applications of EFT 
to nuclear physics and to other systems which exhibit shallow bound states. 

The "pionless" EFT discussed here has its historical roots in the effective range expan- 
sion and in the model-independent approach to three-body physics 0. It considerably 
extends those approaches by making higher order calculations systematic, by including exter- 
nal electroweak currents, and by providing rigorous error estimates. A number of technical 
issues like gauge invariance and the inclusion of relativistic corrections are also much simpler 
in this approach. To gain perspective, it is useful at this point to briefly review what has 
been accomplished towards the "pionless" EFT program. In the two-nucleon sector, the 
modification of the power counting needed to deal with the fine tuning alluded to above and 
the resulting non-perturbative nature of the low-energy expansion is well understood. It has 
been phrased in elegant terms using dimensional regularisation which allows for analytical 
calculations. A large number of processes were analysed, some of them to very high order 
(see Refs. 0, H] and references therein). Typically, most of the low-energy constants 

required are obtained from the plethora of available nucleon-nucleon scattering data. At 
some but usually high order, new two-body operators appear which describe effects due 
to "exchange currents", "quark effects", etc. They contribute typically in the range of a 
few percent to the final result. The coefficients of these operators are not constrained by 
nucleon-nucleon scattering data, but have to be fit to experiments in which the interaction 
of external currents with the two-nucleon system is probed. The low-energy expansion con- 
verges very nicely, the computations are not excessively complicated, and the only obstacle 
in obtaining arbitrary accuracy is the question whether there are enough data available to 
fix the low-energy constants. 
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In contradistinction, we are in the three-body sector still at a more primitive stage. For 
most channels of the three- nucleon system (without external currents), the approach used 
in the two-nucleon sector can be extended in a straightforward way. Analytical calculations, 
however, cannot be performed any more, as the basic Faddeev equation is no longer soluble 
in closed form. All these channels are repulsive, so that consequently no bound states exist, 
and the only observables are the nucleon-deuteron phase shifts which have been computed 
to third order []^, |, ^, |TU|. The only operators not determined by nucleon- nucleon data con- 



tributing to nucleon-deuteron scattering are three-body forces which are highly suppressed 
in these repulsive channels. In contradistinction, the ^He-^H channel exhibits the effect of 
the two-nucleon fine tuning much more dramatically. This leads to a very unusual situ- 
ation. On a diagram-by-diagram basis, the graphs contributing at leading order (LO) to 
nucleon-deuteron scattering are finite. This seems to agree with the expectation, based on 
naive dimensional analysis, that the three-body force, introduced to absorb any divergence, 
does not appear at leading order. However, it turns out that the ultraviolent behaviour 
of the equation describing the three-body system is very different from the behaviour of 
its perturbative expansion. In particular, a non-derivative three-body force is required to 
absorb cutoff dependences already at leading order. This immediately voids the use of the 
effective range expansion. The precise value of the three-body force is not determined by 
nucleon-nucleon data and has to be determined by some additional three-nucleon datum. 
These unusual properties of the Faddeev equation and the existence of a new free parameter 
in the three-body sector was recognised early on |1T[, and the renormalisation theory inter- 



pretation was given in |]T2[. These complications hindered the extension of a consistent and 



economical power counting to the three-body sector, and currently only leading and next- 
to-leading order (NLO) calculations of the doublet S-wave neutron deuteron phase shift and 



of the corresponding triton binding energy are available |T3|, |T^. To attain the precision 
necessary to be of phenomenological interest, one has to go at least one order beyond that. 

This paper is devoted to re-formulating the renormalisation issues in the ^He-^H sector 
in order to extend the power counting to all orders. In the following section, we do so first 
in a bosonic model to avoid unrelated complications coming from the spin-isospin structure. 



In Sect. PH we consider the three-nucleon system, and compute then the neutron-deuteron 
phase shift in this channel at next-to-next-to-next leading order (NNLO) in Sect. There, 
we also compare the results with the available phase shift analysis [TH], and with potential 



model calculations [16|. This also allows us to discuss the convergence of the low-energy 
expansion in order to obtain a credible error estimate. We close by conclusions, outlook and 
an appendix. 



II. BOSONIC INTERLUDE: POWER COUNTING TO ALL ORDERS 

Let us first consider the renormalisation of higher order corrections in the simpler setting 
of the model of three spinless bosons of mass M whose two body system exhibits a shallow 
real bound state. This avoids the technical complications with spin and isospin degrees of 
freedom. The Lagrangean is thus 

^ = (i^o + ^) ^ + t'"^'' + T (^^^^"^ " V) V)^' + H.c.) 

+ ^i^r + ^ ((^n V - V) v)^v' + H.c.) + . . . , 
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where the terms not exphctly shown are further suppressed. A more convenient form of 
writing the same theory is obtained by introducing two dummy fields d and t with the 
quantum numbers of two and three particles (referred to in this section as the "deuteron" 
and "triton") 



^ \ ° 2M J ^ V 4M J V 6M M J 



[d^^ + H.c.) - uj [tU^ + H.c.) + 

v2 



(2) 



The equivalence between the two Lagrangeans (|1]) and (0) is established by performing the 
GauBian integral over the auxiliary fields d and t and disregarding terms with more fields 
and/or derivatives, that are suppressed further p. However, there is a subtlety. In order 
to eliminate time derivatives and express the result of this integration in a manifestly boost 
invariant way as in (|^), one has to perform a field redefinition 

(3) 

The relation between the constants in the two Lagrangeans is then given by 

29' 



Co 
Q 

Do 
D2 



A ' 
^ 4MA2 ' 



(4) 
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We can identify three scales in the problem: Q is a typical external momentum, R is the 
short distance scale in position space beyond which the EFT breaks down, and 7 = \/JAB^ 
the anomalously small binding momentum of the shallow two boson bound state of binding 
energy B2. This scale is, up to effective range corrections, equal to the inverse scattering 
length in the two-body channel. We thus perform an expansion in powers of two small 
dimensionless numbers, QR -C 1 and 7-R <^ 1, while keeping terms with arbitrary powers 
of (5/7 ~ 1. The corresponding power counting has been extensively discussed elsewhere 
and will not be repeated here, see e.g. 0, ^ and references therein. Its details depend on 
the particular regularisation and renormalisation scheme used. In this work, we use a sharp 
momentum cutoff to regulate divergent integrals^. The result of these considerations can be 
summed by stating which diagrams contribute at each order in the low-energy expansion. 

The bare deuteron propagator is given by the constant term i/A. At leading order (LO), 
the d field propagator is given by the sum of all graphs made up from the interaction 
proportional to g in with an arbitrary number of loops, see Fig. |l]. At next-to-leading 
order (NLO), one includes one perturbative insertion of the deuteron kinetic term. At 



At present, it is not known how to extend dimensional regularisation, which is the scheme used in two-body 
calculations, to the three-body sector. 
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FIG. 1: Top: Re-summation of the bare deuteron propagator (thick sohd hne) into the dressed 
deuteron propagator (double hne) at LO by dressing with two-nucleon bubbles. Bottom: The NLO 
and NNLO corrections to the deuteron propagator. The cross denotes an insertion of the deuteron 
kinetic energy operator. 



next-to-next-to- leading order (NNLO), two insertions of the kinetic energy enter. At higher 
orders, terms with more derivatives not explicitly shown in (Q) appear. We thus arrive at 
the deuteron propagator in the low-energy expansion 



(5) 



1 



-7+ J-Mpo + ^-ie 



j^ + Mpo-S^ 



-7 + V-Mpo + ^-ie 



(72 + Mpo - 



-7 + V-Mpo + ^-ie 



+ 0(p3/^3^ 



Air 



where p is the effective range, defined by performing the effective range expansion around 
the bound state. The effective range is assumed to be of natural size, p ~ i?. We also 
defined the rescaled propagator V{pQ,p) for future convenience and quote the two particle 
propagator when the effective range is summed to all orders, but the other corrections in 
the effective range expansion are left out: 



1^{Po,p) 



1 



-7 



-Mpo + S--ie + 



(72 + Mpo - 



(6) 



The two particle scattering amplitude is obtained from the d field propagator by multiplying 
it by g^. After that, we recognise the familiar form of the scattering amplitude in the effective 
range expansion, which had allowed us to determine the low-energy constants Cq, C2 in terms 
of 7 and p in (||/^. 

Similarly, the leading order three-particle amplitude includes all diagrams build up of the 
leading two-body interactions, i.e. the ones proportional to Cq in (P or g and A in (Q). The 
NLO amplitude includes the graphs with one insertion of the deuteron kinetic energy, the 
NNLO amplitude includes those with two deuteron kinetic energy insertions, and at higher 
orders new interactions not explicitly shown in (|ll/0) appear. But that is not all. There also 
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appear graphs containing the three-body interactions proportional to Do, D2 in (|T|) or Q and 
the triton kinetic energy in (^. It has been estabhshed that the leading three body force, 
i.e. the one proportional to Dq in (|T]) or Q in appears at LO and NLO |]12|, |13|, [l^. We 
will re-derive this fact in a slightly modified way which permits us also to find the proper 
generalisation to all orders in the expansion. 

Our argument proceeds as follows. We include a specific three-body force term at a given 
order n of the expansion if and only if that term is needed to cancel cutoff dependences 
in the observables which are stronger than cutoff dependence from the suppressed terms of 
order (Qi?)""*"^. The rationale we follow is thus simply a consequence of the fact that physical 
observables cannot depend on the arbitrary value of the cutoff and, from the point of view 
of the EFT, cutoff dependences coming from loops and from the values of the low-energy 
constants cancel order by order in the expansion. 

Before we implement this idea, let us make a few comments about the economy of com- 
puting diagrams at high orders. There is an infinite number of diagrams contributing at 
leading order, as shown in Fig. ^. We sum them up by solving the Faddeev integral equation 
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FIG. 2: Re-summation of the naive LO three-body equation, without three-body interactions (top), 
into the corresponding Faddeev integral equation (mid) for the Boson case. Bottom: The same 
including three-body interactions with strength Hq{A), denoted by the dot. Remaining notation 
as in Fig. |l|. 



for point-like interactions. In the case where only two-body interactions are present 

This re-summation is necessary since 



this 

according to the 



equation was first derived in |T7 
power counting - all these diagrams contribute equally to the final amplitude. The NLO 
contribution can then be obtained by perturbing around the LO solution with the deuteron 
kinetic energy operator, as shown diagrammatically in Fig. ^. This approach was taken be- 
fore and is easily implemented However, it suffers a drawback at higher orders because 
one must then insert the kinetic energy operator two or more times. This is very cumber- 
some to do numerically, particularly because of the need to compute the full off-sheU LO 
amplitude before inserting the deuteron kinetic energy terms, see Fig. |^. A more practical 



6 



FIG. 3: The graphs to be computed when inserting the higher order corrections to the deuteron 
propagator perturbatively. Top: NLO. Bottom: NNLO. The bottom right diagram makes the 
knowledge of the full off-shell LO amplitude mandatory. Notation as in Fig. |^. 



method would be to re-sum all range corrections into the deuteron propagator and solve a 
modified integral equation. This approach has its own drawbacks, including the existence 
of spurious bound states in the re-summed deuteron propagator, and was shown to fail 
miserably for the phase shifts 



1^ 



We choose a middle ground: First, expand the kernel of the integral equation perturba- 
tively, and then iterate it by inserting it into the integral equation. That means, we perform 
a partial re-summation of the range effects: We use the deuteron propagator in (||) in the 
integral equation, which is then solved numerically. Figure ^ shows a graphical representa- 
tion of the NLO version. This procedures includes some graphs of higher order, for instance, 
the NLO calculation includes diagrams like the bottom right one in Fig. |^, but not diagrams 
like the bottom left one in the same Figure. The partial re-summation arbitrarily includes 
higher order graphs and consequently does not improve the precision of the calculation. But 
the additional diagrams are small in a well-behaved expansion, so the precision is not com- 
promised, either. On the other hand, this prescription is much easier to handle numerically 
than the two other methods. In a nutshell, the only necessary re-summation is the one 
present at LO. The partial re-summation of range effects is made only for convenience. 
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FIG. 4: The integral equation for the boson case when the kernel is expanded to NLO. Notation 
as in Figs. || and 
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The integral equation describing particle-deuteron scattering one obtains by including the 
re-summation discussed above is hence a minor modification of the one derived in [12, l7 . 
Its NLO version is graphically represented in Fig. ^. The integral equation for the half 
off-shell amplitude up to and including NNLO reads: 

A ^ 

tA(p, k) = /C(p, k) + Hip, A) + ^Jdq q'V{E - j^.q) [/C(p, q) + n{E, A)] t^{q, k) , (7) 



where 



n{E,A) = 7^o(A)+^2(^,A) , (8) 
= ^ (ME + ,^) = (ME + y) , (10) 



^(P,.)^^ln( ^^^^!^!!^ ) , (11) 
pq — pq + q^ — ME 



and E = 3/c^/(4M) — is the total centre-of-mass energy. The on-shell scattering 

amplitude T{k) is given by Zt\{k^k) , where the deuteron wave function renormalisation 
factor is given by 

z=^{i+ip+iipr+...) . (12) 

and is related to the phase shift 6 in this channel by the usual formula 

Tik)^Zt.(k,k)^-^-^^. (13) 

As (|^) contains the cutoff A, the solution tx{q) for arbitrary Hq, H2 depends in general on 
the cutoff, too. Since the low-energy amplitude can not depend on the cutoff, any exphcit 
dependence on the cutoff has to be cancelled by the cutoff dependence implicit in Hq, H2, 
etc. We determine the cutoff dependence of the three-body forces by imposing this condition 
order by order in 1/A. In this process, we also determine which three-body force appears at 
every order of the expansion. 

Let us first consider the effect of changing the cutoff from A to A' ~ A. The equation 
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satisfied by the amplitude ^a' analogous to can be written as follows: 



J dq {6{p - q) q'V{E - q)IC{p, g)) = /C^ + HiE^ 



^l/Q2 ^l/Q2 



A 

+ ^ J dq q'V{q)n{E,A)t^,{q) 



A/(QA) 



+ ^ J dq q'V{E - ^,q){n{E,A') -niE,A))Mq) 



+ n{E,A')-n{E,A) 

~1/A2 ~1/A2 



A' 

^2 



+ ^ I dq q'V{E - ^,q) (^}C{p,q) +niE,A'))tj,,{q) (14) 

""^^i/IqaT ~i/[qa) 

The point of this rewriting is to show explicitly as in the last three lines the difference 
between the equations satisfied by t^ip) and t\/{p). Now we take the cutoffs A and A' to 
be below the short distance scale (5^A<^l/r~ 1/R, but higher than the infrared scales 
Q J k ^ p. The asymptotic behaviour of t{q) is known to be of the form t{p ^ Q) ~ 
sin(solnp/A)/(Qp) ]19[ (see Appendix), so t{p ~ Q) ~ l/Q^ and t{p ~ A) ~ 1/{QA). As 



one knows the asymptotic behaviour of tx{p) and assumes Hq ^ H2 ^ 1, one can estimate 
the size of each term in (ffl) as shown in under-braces. Here, loop momenta q count as Q 
in ultraviolet finite integrals but count as A in divergent integrals, since the integrals are 
dominated by momenta of the order Q and A, respectively. 

The estimate suggests that the terms on the right hand side describing the contribution 
from loop momenta between A and A' as well as the contribution from the three-body forces 
are small compared to the leading ones on the left hand side, the former being suppressed by 
at least one power of Q/A. Therefore, one naively expects them to be negligible at leading 
order. If one indeed discards them, the resulting equation for t^'ip) becomes identical to (^ 
and one has to conclude that t\/{p) ~ t\{p) + 0{Q/ A). 

The argument above is, however, incorrect. This is due to the fact that left hand side of 
(|l^) is singular in the limit A — 00. More precisely, in this limit, there is a zero eigenvalue 
for the operator on the left hand side of ( |14|) acting on t^(p): 

00 

j dg [5{p -q)- ^q'V{E - q)IC{p, q)^ a{p) = . (15) 



We refer the reader to Ref. [11 1 which proves this statement rigorously and also makes it 



clear that the spectrum of the operator in ([15D is continuous. One important point to notice 



is that the existence of the zero mode is an ultraviolet phenomenon, independent of the 
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infrared scales like E or 7. In the simpler case 7 = 0, the eigenfunction corresponding to an 
eigenvalue A can be expressed analytically at LO as 



a.(p) = si„^-i.l„^yj^+yi-^jj (16) 

with eigenvalue 

8 sir^ (17) 

v/3s cos(7rs/2) ^ ^ 

and —1 < Re(s) < 1 The zero mode ao{p) corresponds to s(A = 0) = iso = i 1.00624. . .. 
The existence of the zero mode shows that for A — 00, the operator on the left hand side of 
(^) has no inverse. Therefore, the solution to the integral equation is not unique since one 
can always add ao{p) to a given solution and obtain another valid solution. 



For finite A, general theorems guarantee that the spectrum of the operator is discrete pO 
Eq. (^ then defines the amplitude uniquely. However, as A is increased the eigenvalues get 
closer to each other and, for a generic value of the cutoff, the eigenvalue closest to zero is of 
the order 0{Q/A). Since the operator in the left hand side of (|^) is nearly singular, terms 
of order 0{Q/A) on the right hand side of (|1^ can have an effect on t/\^{p) of order 0{1). 
This "hypersensitivity" to cutoff dependences is at the heart of the apparently paradoxical 
fact that the solution of (0) with H = is cutoff dependent, even though the contribution 
of the high momentum modes are suppressed by 1/A (that is, are ultraviolet finite). In 
general, for the amplitude to be cutoff independent up to terms of order (9((Q/A)") , the 
right hand side of ( [T^ ) has to be cutoff independent up to terms of order (^((Q/A)""'"^) , 
i.e. one order higher than in the standard cases. Notice that the unexpected dependence 
of the low-energy amplitude on high momenta is due to two factors: the harder asymptotic 
behaviour of t{p) ~ 1/p compared to the perturbation theory expectation t{p) ~ 1/p^, and 
the fact that the leading terms in (p!4D form a nearly singular equation that amplifies any 
remaining term by one power of A/Q. 

We now determine the three-body forces by demanding order by order in Q/A 

A ^ 

n{E,A) + ^ j dq q'ViE - ^,q)[}Ciq,p) +niE,A)]tAiq) = const , (18) 

i.e. the amplitude must be cutoff independent. This (possibly energy dependent) constant 
can without loss of generality be set to zero by absorbing it into a re-defined three body 
force 

n{E, A) ^ n{E, A) . (19) 

,2 



l + ljdq q'V{E-^,q) t^{q) 



Imposing (|T^) with const. = order by order in the low energy expansion, we can determine 
which three-body forces are required at any given order, and how they depend on the cutoff. 
As we are interested only in the UV behaviour, i.e. in the asymptotics, and no IR divergences 
occur, we can safely neglect the IR limit of the integral. 
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At leading order, 7 <^ A, and we can approximate 



^ ^ + • • • 



(20) 



t(g)^A(A:/A,7/A)^^"^^°^"^^/^^U... . 

where a normalisation factor A was included into the asymptotic form of t{q). The essential 
observation is that the terms of order 0{Q/A) are independent of E (or k) and hence can be 
made to vanish by adjusting only Hq and not H2 or any of the higher derivative three-body 
forces. In fact, using the substitutions of ( pOD in (0) and neglecting higher order corrections 
for naturally-sized Hq, we find 

_ sin^lnf + arctango] ^^^^ 
sin[so In ^ — arctan sq] 

as obtained previously [|T^]. 

At NLO, we impose (|I^) to order (9(1/A^). Thus, we expand one order higher than (^) 

/C(g,p)~4+- 

The sub-leading asymptotics at NLO depends on whether A<^l/rorA^l/r. In the 
first case, there is only a perturbative change from the LO asymptotics, and we can find it 
analytically (see Appendix) 

t{q) ~ A{k/A) Pi^[^oln(g/A)] ^ i^^i ^ gj^j^^ Hq/A) + arg(Si)] 

^ - (23) 

sm[soln(g/A) + arg(5_i)] 

+ -D-l 07 h . . 

g2 

Again we see that, to this order, no terms dependent on E or k are present. Only the Ho{A) 
term is necessary for renormalisation. Plugging (p^/|23|) back in (|18|) and keeping only the 
terms of order l/A^, one can derive the sub-leading correction to (pi]), but the expression is 
long and of no immediate use. 

A new three-body force enters at NNLO. We require now that the terms of order (Q/A)^ 
on the right hand side of (|I^) should vanish, which guarantees that the amplitude is cutoff 
independent up to terms of order (Q/A)^. At this order though, there are terms proportional 
to ME arising from expanding the kernel in powers of q/A and V ME/ A. Indeed, such 
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corrections appear only at even orders. Only a three-body force term which also contains 
a dependence on the external momenta k"^ and ME (both equivalent on-shell) can absorb 
them. That means that H2 appears at NNLO. In addition, there is a large number of 
momentum independent terms that renormalise Hq and provide the NNLO corrections to 
(prp. The same pattern repeats at higher orders. Every two orders, a new three-body force 
term appears with two extra derivatives, guaranteeing that the new three-body forces are 
analytic in ME oc on shell. 

If we choose A > 1/r, the argument is changed only by the change in the asymptotic 
form of t{p). This affects the explicit forms for the running of Hq, H2, but not the fact that 
they are necessary and sufficient to perform the renormalisation. 



III. THREE NUCLEONS IN THE TRITON CHANNEL 

It is not difficult to include the spin and isospin degrees of freedom and generalise the 
discussion in the previous section to the case of the three-nucleon system. As discussed 
elsewhere, the ^Si channel - to which ^He and belong - is qualitatively different from 
the other three-nucleon channels. This difference can be traced back to the effect of the 
exclusion principle and the angular momentum repulsion barrier. In all the other channels, 
it is either the Pauli principle or an angular momentum barrier (or both) which forbids the 
three particles to occupy the same point in space. As a consequence, the kernel describing 
the interaction among the three nucleons is, unlike in the bosonic case, repulsive in these 
channels. The zero mode of the bosonic case, discussed above, describes a bound state since 
it is a solution of the homogeneous version of the Faddeev equation. As such, it is not 
expected to appear in the case of repulsive kernels and, in fact, it does not. For the '^He-^H 
channel however, the kernel is attractive and, as we will see below, closely related to the one 
in the bosonic case. Consequently, the arguments in the previous section are only slightly 
modified for the case of three nucleons in the triton channel. 

The three-nucleon Lagrangean analogous to (^ is given by 

+ ('^° + ^ + M + ^) ^ " (rff (AT^P^AT) + H.c.) - gt (f^N^P^N) + H.c. 
- Us {t\T^N)dj + H.c.) - ut {t^{a'N)di + H.c.) + . . . , 

(24) 

where A^ is the nucleon iso-doublet and the auxiliary fields t, df and dl carry the quan- 
tum numbers of the ^He-'^H spin and isospin doublet, ^5*0 di-nucleon and the deuteron, 
respectively. The projectors and P"^ are defined by 



1 (25) 



where A = 1,2,3 and i = 1,2,3 are iso-triplet and vector indices and (a*) are isospin 
(spin) Pauli matrices. 
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When the auxihary fields are integrated out of (p4D, three different kinds of three-body 
force terms without derivatives seem to be generated: 

^3body = {N^P'Ny{N^a'a'N){N^P'N) 

= -{N^P^Ny{NW^a'N){N^P'N) ^ ^ 

oc {N^Nf . 

This last form for the three-body force is explicitly Wigner-SU (4)-symmetric, i.e. it is sym- 
metric by a generic unitary transformation of the four fields corresponding to the 2 (spin) 
x2(isospin)= 4 nucleon states E^ and corresponds to the singlet and triplet effective range 



parameters to be identical. As observed before, the term in ( pB]) is the only possible three- 
body force without derivatives. Thus, the different couplings u!s,uJt are redundant and we 
can take them to be determined by only one free parameter 

A2 ~ Mg^ n ~ Mg^ n ~ Mg^gt n ' ^ ^ 

Contrary to the terms without derivatives, there are different, inequivalent three-body force 
terms with two derivatives. As we will see below, however, only the S'[/(4)-symmetric 
combination is enhanced compared to the naive dimensional analysis estimate and appears 
at the order we work. We keep then only this special combination. This amounts to choosing 

A4 ~ M2^2^2 - ]vf2^2^2 - M'^g.gt^'^ ' ^ ' 

With these choices, we have six parameters left at NNLO: IS.s,gs (determined by the scat- 
tering length and effective range in the ^5*0 two-nucleon channel), At,gt (determined by the 
deuteron pole and effective range in the deuteron channel), Hq and H2 (to be determined 
by three-nucleon data). 

The derivation of the integral equation describing neutron-deuteron scattering has been 



discussed before |T3|, ^7j. We present here only the result, including the new term generated 
by the two- derivative three-body force. Two amplitudes get mixed: tg describes the dt+N — *■ 
ds + N process, and describes the dt + N~^dt + N process: 

A 

Up, k) = \ [3/C(p, k) + 2H{E, A)] + ^ y dg [P,(g) [/C(p, q) + 2n{E, A)] t,(g) 



+A(g) [^lC{p,q) + 2H{E,K)]tt{q)] 

A 

tt{p, k) = \ [/C(p, k) + 2n{E, A)] + ^ y dg q^ [A(g) [/C(p, q) + 2n{E, A)] tM) 



+V,{q) [3/C(p, q) + 2H{E, A)] U{q)] (?9) 

where Vs,t{q) = 'T^s,t{E — j^,q) are the propagators analogous to (|]). For the spin-triplet 
S-wave channel, one replaces the two boson binding momentum 7 and effective range p by 
the deuteron binding momentum 7^ = 45.7025 MeV and effective range pt = 1.764 fm. 
Because there is no real bound state in the spin singlet channel of the two-nucleon system, 
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its free parameters are better determined by the scattering length = l/^^s = —23.714 fm 
and the effective range = 2.73 fm at zero momentum, i.e. instead of 

VM = , ^ + (30) 



(-7. + \/l{p'-k^)+l?) ^ ^ {-Is + ^JW-k^)+li) 

The neutron-deuteron J = 1/2 phase shifts 5 is determined by the on-shell amphtude 
tt{k, k), muhiphed with the wave function renormahsation for the deuteron as given in (|T2p ^ 

TW^Z,(M)^|^^^. (31) 

The renormahsation of ( p9D is better understood after introducing the variables t+ = tg + tt, 
t_=ts- tt, V+ = (V, + Vt)/2 and V_ = (P, - Vt)/2, in terms of which (M) reads 



t+{p,k) = }C{p,k)+n{E,A) 

A 

+ ^ J dq q' []Cip,q)+n{E,A)] [I)+(g)t+(g) + 2)_(g)t_(g)] (32) 



A 

t4p,k) = i }C{p,k) -IJdq q'lC{p,q) PMUq) +T^-{q)t+{q)] • (33) 



The essential observation is that in the Wigner-5'f/(4)-limit (P_ = 0), the two equations in 
(p9|) decouple and that, in the ultraviolet regime where the differences between 7^ and jt 
disappear, the SU{4) limit is recovered. The amplitude t+ satisfies the same equation (□) 
as the boson case above, and t_ satisfies the same equation as the '^Ss amplitude [|13], \l% 
which does not require any renormahsation at this order. Thus, the arguments given for the 
renormahsation in the bosonic case applies with only minor changes for the ^Si channel. 
Let us now look at these changes in detail. First, we note that the asymptotic behaviour 

2 

of t-(p, k) is driven by the term containing V^{E — g)t+(g): 



00 

t^{p,k) ~ y"dg q'^}C{p,q) 



7s - 7t sin[soln(g/A)] 



(34) 

sin[so ln(p/ A) + arg B_i] 
~ (7s - It) 2 ■ 

We can separate (^) and (|33| ) as in ([T^), retaining on the left hand side the piece proportional 
to P+. Eq. ( p^ ) has the same zero mode we found in the bosonic case so, just as in the 
bosonic case, the amplitude is sensitive to the pieces of order 1/{QK) on the right hand 
side already at leading order. On the other hand, thanks to a crucial change of sign in the 



The mixing with the ^Di channel vanishes at the order we consider here. 



14 



kernel, (|33D does not have a zero mode and does not present this enhanced sensitivity on 
1/A suppressed terms. The terms of order 1/A or higher can be estimated by approximating 
the kernel in the ultraviolet as in (pO|/p^^ 



V_{E 



2M'^^ 2 



4 1 4 7^ + 7t ps + Pt 
3 o2 + 2 



2M 



47 
3 



It 



+ 



Ps- Pt 



(35) 



By combining the asymptotic behaviour of t+, t_ and (^), one can find the contribution of 
the high momentum modes (p ~ A) from the different terms in (^/ |33|) . Starting from (p2|), 
we see that first, V^t^ is of order 1/(QA); second, it is identical to its bosonic counterpart; 
and third, it is absorbed by the same Hq as in (^). The same term also generates a term 
~ 1/A^ (equal to the one in the bosonic case after the substitution 7 — (7^ + 7t)/2,r — >■ 
{Ps + Pt)/2 ), and one scaling like ~ ME / {Qh?), absorbed into Hq and respectively, 
in close analogy to the bosonic case. The term contributes with a term of order 

^ {is ~ 7t)^/(QA'^), which is absorbed into Hq. As mentioned before, (^) does not exhibit 
the enhanced sensitivity to terms suppressed by powers of 1/A. Thus, to guarantee that 
t_ is cutoff independent up to terms of order 1/((5A"), it is enough to make (|33[) cutoff 
independent up to terms of the same order. Both terms "D+t- and contribute with 

terms starting at order (7^ — 7t)/((5A^). However, it is well-known that the only point-like 
three-body force without derivatives is necessarily SU (4)-symmetric. In the renormalisation 
group sense, the hence irrelevant (i.e. spurious) interaction induced at finite A is irrelevant 
and serves only to repair violations of Wigner's S'[/(4) symmetry by our choosing a cutoff, 
which is equivalent to smearing out the three-body operator. In closing, we stress that the 
approximations in ( |55| ) are performed only in the ultraviolet regime, when analysing the 
contribution of order ~ 1/(QA) to the integral equation. This analysis determines only 
which three-body forces are needed to render cutoff invariant results. The full equation with 
7s 7^ It, Ps 7^ Pt has to be solved to compute the phase shifts, and to adjust the three-body 
forces. 



IV. NUMERICAL RESULTS FOR ^Si NEUTRON-DEUTERON SCATTERING 



We numerically solved the coupled equations in (|2^) up to NNLO using besides the 
values given above the iso-spin averaged nucleon mass M = 938.9 MeV. The com- 
putational effort in solving Faddeev equations for separable potentials, as it is the case 
here, is trivial, and the reader is invited to download the Mathematica codes used from 
[http : / /www-nsdth . Ibl . gov . To deal with the deuteron pole and the logarithmic singu- 



larity (which is integrable but generates numerical instabilities), we follow Hetherington 
and Schick to first solve the integral equation on a contour in the complex plane and 
then use the equation again to find the amplitude on the real axis. We determine the two- 
nucleon parameters from the deuteron binding energy and triplet effective range (defined by 
an expansion around the deuteron pole, not at zero momentum) and the singlet scattering 



^ For the sake of the simphcity of the argument, we take here ps — which is in general true up to higher 
order corrections, in order to use the same effective range expansion of singlet and triplet NN . 
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length and effective range (defined by expanding at zero momentum). Finally, we fix the 
three-body parameters as follows: Because we defined H2 such that it does not contribute 
at zero momentum scattering, one can first determine Hq from the ^Si scattering length 



03 = (0.65 ± 0.04) fm [23|. At LO and NLO, this is the only three-body force entering, but 
at NNLO, where we saw that H2 is required, it is determined by the triton binding energy 
i?3 = 8.48 MeV. The cutoff is varied as discussed below. 

The error in the determination of the ^Si scattering length of about 7%, last measured 
1971, dominates our errors at NNLO. A better experimental determination of this parameter 
is not only a necessary condition for the improvement of the accuracy of our method. It is 
also important information for fixing any three-body force in sophisticated models of nuclear 
physics, given that there is only scarce information on phase shifts in the triton channel. 




FIG. 5: Logarithm of the cutoff variation of phase shifts between A = 200 MeV and A = 600 
MeV as function of the centre-of-mass momentum. The short dashed, long dashed and solid line 
correspond respectively to LO, NLO and NNLO. The dots correspond to a "NNLO" calculation 
without the three-body force H2. 

The importance of isospin- violating effects can be estimated to be of higher order. In fact, 
the splitting between the ^5*0 scattering lengths between neutron-neutron («„«) and neutron- 
proton ( be estimated as (l/a„„, — l/a„p)/Q, where Q is the typical momentum in 
the system. Taking Q ~ ^fjABl this is a correction of about 2%, smaller than the cutoff 
variation we see at NNLO. For future high precision calculations, however, their inclusion 
is certainly required. 

Unless otherwise stated, all calculations were performed with the cutoff at 600 MeV. 
In the phase shift plots, we include the range of the NNLO result as the cutoff is varied 
between 200 and 600 MeV. Albeit it neglects other sources of error at low cutoffs, this is a 
reasonable estimate of the errors, especially at high momentum. A definite answer for the 
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size of the errors, however, can only be given by a higher-order calculation. The lower limit 
of this variation was chosen to be of the order of the pion mass. The upper limit is rather 
arbitrary since beyond A ^ 400 MeV, the phase shifts are essentially cutoff independent. 
The phase shifts are cutoff independent up to the order the calculation is valid. This is 
demonstrated in the plot in Fig. |^, where we show the relative difference in k cot S caused by 
varying the cutoff between A = 200 MeV and A = 600 MeV at different orders. Two points 
are worth noticing. First, the cutoff variation decreases steadily as we increase the order of 
the calculation and is of the order of (fc/A)", (7/A)", where n is the order of the calculation 
and A = 200 MeV is the smallest cutoff used. Second, the errors increase with increasing 
momentum, as one would expect from the fact that there are errors of order (fc/A)". The 
slope does not change significantly as we go to higher orders due to the (7/A)"' terms. We 
also show in the same figure the same cutoff dependence of phase shifts but now calculated 
at NNLO without the three-body force H2. It shows that not much improvement is obtained 
over the cutoff dependence of the NLO calculation, lending credence to the power counting 
discussed above. 

The phase shifts at LO, NLO and NNLO are shown in Fig. together with the result 
of a phase shift analysis [|15| and results from numerical calculations using the two-nucleon 



Argonne V18 potential and the Urbana IX three-body force [16| . The phase shifts converge 



steadily in the whole region of validity of our theory, up to centre-of-mass momenta of the 
order of 120 MeV. At a momentum of the order of y/MB^ ~ 90 MeV (-B3 the triton binding 
energy) the cutoff variation is of order of 24%, 12% and 4% respectively at LO, NLO and 
NNLO. Convergence by itself and to the results of the more sophisticated potential model 
calculation is achieved, albeit less pronounced in the imaginary part of the phase shift. 

The comparison of our results with those obtained through potential models must however 
be done with care. Results obtained by the use of effective field theory have an universal 
validity in the following sense. Any model with the same experimental values for the low- 
energy constants as the ones we used to fix the low-energy constants of the effective theory, 
shares also the same prediction for other observables, up to the accuracy of the order of the 
EFT calculation. In our case, we expect thus that a model like the Argonne V18 two-body 
force complemented by the Urbana IX three-body force which predicts or fits the triton 
binding energy and low-energy two-body phase shifts, should give phase shifts agreeing 
with the effective theory result at the 4%, and this is indeed found. The agreement may 
disappear at higher order: The effective theory result will depend on low-energy constants 
determined by three-body observables that do not enter in the potential model determination 
and, through this information, can acquire higher accuracy. This is especially true when 
electroweak currents are considered. 

A next-order calculation with an accuracy of about 1% is straightforward and must 
include the mixing between the 1 and ''Di partial waves, the contribution of two-body p- 
waves and isospin-breaking effects. However, although it does not involve any new unknown 
low-energy constants, it will not provide improved precision, since of the parameters used 
to fix our constants at the lower orders, the scattering length in the doublet S wave channel 
is not known with sufficient accuracy and would dominate the errors. 

The universality of the EFT predictions mentioned above is well illustrated by the Phillips 
line It is a well-known empirical fact that different models having the same low- 

energy two-body physics may have widely different three-body physics. The variation on the 



We thank A. Kievsky for providing us with his results for the phase shifts at our request. 
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k [MeV] 



FIG. 6: The neutron-deuteron ^Si phase shift at LO (dotted Une), NLO (dashed Une) and NNLO 

2 

(dark band). Top: Real part. Bottom: Imaginary part. The band in the NNLO curve corresponds 
to the change caused by a variation of the cutoff from A = 200 MeV to A = 600 MeV. The dots are 
from the only available direct phase shift analysis of the triton channel [p!^ ], which dates from 1967, 
and the crosses are the results obtained with the Argonne V18+Urbana IX two and three-body 
forces [16|. 
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FIG. 7: The Phillips line predicted at NLO. The dots correspond to the predictions for the triton 
biding energy and doublet scattering length in different models with the same two-body scattering 
lengths and effective ranges as inputs [^] . The dotted and dashed lines are the EFT results at LO 
and NLO, respectively. The cross is the experimental result. 

three-body physics is, however, correlated: One single three-body observable, say the triton 
binding energy, determines all the others, up to small corrections. In the EFT language, 
this is a reflection of the existence of the parameter Do (or Ho) already at leading order 
which is not determined by two-body physics but whose precise value at fixed cutoff must 
be determined from another three-body datum. Neither Dq nor A has an immediate physical 
significance, but once it is determined at a given cutoff, its cutoff dependence is fixed by 
(p!8|). Therefore, as either of these parameters is varied as the other is kept fixed, the results 
of different potential models are reproduced. The basic physics leading to this effect was 
explained by Efimov in ^ . We here exemplify it by showing a plot of the predictions for the 
triton binding energy versus the doublet scattering length in a variety of models in Fig. 
together with the EFT result at LO and NLO. The NLO result includes the (convenient but 
unnecessary) re-summation of higher order terms discussed above. We note that our NLO 
result is very similar to Efimov's that was obtained by performing perturbation theory on 



the effective range 
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V. CONCLUSIONS AND OUTLOOK 



We presented a simple scheme to classify at which order a given three-body force must 
be considered in an EFT approach to the few-body system at very low energies. These 
analytical considerations rest on the premise that a three-body force is included if and only 
if it is necessary to cancel cutoff dependences in the observables at a given order. 

We compared the predictions of this scheme with its numerical implementation in the 
case of a NNLO calculation of the nd scattering phase shifts in the ^Si (triton) channel. 
We confirmed the well-known result that to LO and NLO, only one momentum dependent 
three-body force is necessary for cutoff independence, whose strength can be determined by 
the three-body scattering length. We found that one and only one new parameter enters 
at NNLO, namely the Wigner-S'[/(4)-symmetric three-body force with two derivatives. The 
one additional datum needed to fix it is the binding energy of the triton. Indeed, only the 
S'f/(4)-symmetric three-body forces are systematically enhanced over their contributions 
found from a naive dimensional estimate, with a 2n-derivative three-body force entering at 
the 2nth order. On the way, we also propose a new, computationally convenient and simple 
scheme to perform higher order calculations in the three-body system by iterating a kernel to 
all orders, which has been expanded to the desired order of accuracy in the power counting. 

The phase shifts thus calculated are in good agreement with a partial wave analysis and 
more sophisticated model computations. The result shows convergence. 

Work under way includes the NNNLO calculation of the triton channel, its wave functions, 
the extension to ^He and the pd channel, and the inclusion of electro-weak interactions. This 
will allow comparison to a plethora of data and lead to predictions in the range and to an 
accuracy relevant for astro-physical problems. Our scheme is clearly extendible to any EFT 
in which at least a partial re-summation of graphs at LO is necessary. More formal aspects 
of our method must also be addressed. 
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APPENDIX A: ASYMPTOTIC AMPLITUDE 

In this appendix, we derive the asymptotic form of the half off-shell amplitudes quoted in 
the main text. Consider first the bosonic integral equation in (0) in the regime 7, \/ME <^ 
p -C A. For these values of p, the integral is dominated by values of q of the order of p. We 
can then disregard the infrared scales 7 and V ME as well as the term Hq/A"^ (assuming 
Hq ^ 1). The integral can be extended to infinity by adding a sub- leading piece. First 
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dropping the inhomogeneous term, we find 



t(p) = 



oo 




t{q) . 



A Mellin transformation determines tlie solution to be a power law t{p) ~ p' 
determined by the equation 

Ks) = 1 



1-1 



(Al) 



where s is 



with 



oo 



I(s) = 



V3: 



Ivr 



da; In 



x"^ + X + 1 

— X + 1 



8 sin(f ) 

V3s cos(f ; 



(A2) 



(A3) 



The value of s satisfying ( [A^ ) with the largest real part dominates the asymptotics of t{p). 
These values are the imaginary roots s = ±iso ~ ±i 1.0062 .... Because the inhomogeneous 
term in (|^ behaves like ~ 1/p^ in the asymptotic regime, it is consistent to ignore it 
compared to the asymptotics p^'^o-i Qf i^j^g homogeneous term. The overall normalisation 
of t{p) can however not be determined by this analysis. What is perhaps more surprising 
is that the phase of the solution is not determined, either, since both p^^o-i g^^^j ^-iso-i g^j^-g 
equally acceptable solutions. Sub-leading corrections to the asymptotics can be obtained 
in a similar fashion. For the corrections suppressed by one power of either 7/p or pp for 
instance, we insert t{p) = + 87i?_ip^'*"~^ + Bipp^'^^ in the equation 



t{p) 



dg 



P 



2 1 27 r 
7!g + 3^"^2 



In 



p^ + pq + q^ 
p"^ — pq + q^ 



(A4) 



and keep terms linear in 7 and r. Matching order by order, we find 

a/3 /(iso + n) 



Br, 



4 1 - /(iso + n) 



(A5) 



Taking real combinations of the solutions corresponding to Sq and —Sq, we finally arrive at 
the result quoted in (P5|). Proceeding in a similar way, one can obtain this kind of correction 
to all orders, if necessary. We were however unable to find the sum of this series. 
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